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ABSTRACT 

The astrophysical origin of the r-process nuclei is still unknown. Even the most promising 
scenario, the neutrino-driven winds from a nascent neutron star, encounters severe difficulties 
in obtaining requisite entropy and short dynamic timescale for the r-process. In this study, the 
effect of anisotropy in neutrino emission from a proto-neutron star surface is examined with semi- 
analytic neutrino-driven wind models. The increase of neutrino number density in the wind owing 
to the anisotropy is modeled schematically by enhancing the effective neutrino luminosity. It is 
shown that the neutrino heating rate from neutrino-antineutrino pair annihilation into electron- 
positron pairs can significantly increase owing to the anisotropy and play a dominant role for 
the heating of wind material. A factor of five increase in the effective neutrino luminosity results 
in 50% higher entropy and a factor of ten shorter dynamic timescale owing to this enhanced 
neutrino heating. The nucleosynthesis calculations show that this change is enough for the robust 
r-process, producing the third abundance peak (A = 195) and beyond. Future multi-dimensional 
studies with accurate neutrino transport will be needed if such anisotropy relevant for the current 
scenario (more than a factor of a few) is realized during the wind phase (<~ 1 — 10 s). 

Subject headings: nuclear reactions, nucleosynthesis, abundances — stars: abundances — stars: neutron 
— supernovae: general 



1. Introduction 

The astrophysical site of the rapid-neutron- 
capture nucleosynthesis (r-process), which ac- 
counts for about half of nuclei heavier than iron, 
has been a long-standing mystery. During the last 
decade, the neutrino-heated ejecta from a nascent 
neutron star (neutrino-driven winds, Woosley et 
al. 1994) has been considered to be the most 
promising astrophysical site for the r-process. 
Previous studies show, however, severe problems 
in obtaining requisite high entropy and short dy- 
namic timescale for the production of heavy re- 
process nuclei (Qian & Woosley 1996; Otsuki et 
al. 2000; Sumiyoshi et al. 2000; Wanajo et al. 
2001; Thompson et al. 2001). The general rela- 
tivistic effect for a very compact proto-neutron 
star (e.g., the mass of 2.0 Mq with the radius of 
10 km, Otsuki et al. 2000; Wanajo et al. 2001) or 
a magnetar-like magnetic field strength (Thomp- 



son 2003; Suzuki & Nagataki 2005) have been in- 
voked to increase entropy and reduce the dynamic 
timescale of the winds. It is questionable, how- 
ever, if such physical conditions can be the general 
requirements for the r-process nucleosynthesis. 

Qian & Woosley (1996) have suggested that 
an additional energy input to the neutrino-driven 
wind at between 1.5 and 3 times the neutron star 
radius is efficient to increase entropy and reduce 
dynamic timescale of the wind material. In this 
Letter, it is shown that strong anisotropy in neu- 
trino emission from the proto-neutron star, if it 
exists, acts as this extra energy source and helps 
the r-process. The neutrino-driven wind model 
with spherically symmetric, steady outflow ap- 
proximation is used to obtain the wind trajectories 
(§ 2). A sudden increase of neutrino number den- 
sity in winds owing to anisotropic neutrino emis- 
sion is modeled by enhancing the neutrino lumi- 
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nosity. Nucleosynthesis calculations with the ob- 
tained thermodynamic trajectories are performed 
to demonstrate this effect (§ 3). Finally, a possi- 
ble origin of this anisotropy in neutrino emission 
and some implications of this study are discussed 
(§ 4). 

2. Wind Models with Anisotropic Neu- 
trino Emission 

The wind trajectories in this study are obtained 
using the semi-analytic, general relativistic model 
of neutrino-driven winds explored in Otsuki et 
al. (2000) and Wanajo et al. (2001, 2002). In 
this model, the system is treated as time station- 
ary and spherically symmetric, and the radius of 
the neutron star is assumed to be the same as 
that of the neutrino sphere. The heating source 
that drives matter from the neutrino sphere is 
due to neutrino interactions. Heating is due to 
v e and v e capture on free nuclcons (q V N), neu- 
trino scattering by electrons and positrons (que), 
and neutrino-antineutrino pair annihilation into 
electron-positron pairs (q vv ). Cooling is due to 
electron and positron capture on free nuclcons 
(g e Af) and electron-positron pair annihilation into 
neutrino-antineutrino pairs (q ee )- The rms aver- 
age neutrino energies are taken to be 10, 20, and 
30 MeV, for electron, anti-electron, and the other 
flavors of neutrinos, respectively. The mass ejec- 
tion rate at the neutrino sphere M is determined 
so that the wind becomes supersonic through the 
sonic point. 

The mass and radius of the neutron star are 
taken to be M = 1.4 M & and R = 10 km, re- 
spectively. The neutrino luminosity of one spe- 
cific flavor is assumed to be the same for all other 
flavors, which is taken to be a constant value 
L v = 1 x 10 51 ergss _1 . Note that the assumption 
of a constant L„ is reasonable, since the crossing 
time of a wind over the heating region (< 30 km, 
see Fig. 2) is short enough, ~ 0.1s, compared to 
the decay timescale of L u (a few seconds, e.g., 
Woosley et al. 1994). As explored in previous 
studies, this typical choice of parameter set (with 
isotropic neutrino emission) results in insufficient 
physical conditions (i.e., low entropy and long dy- 
namic timescale) for the production of heavy r- 
process nuclei (e.g., Wanajo et al. 2001). 

In this study, anisotropy in neutrino emission 



is modeled schematically as follows. Given there 
is substantially higher neutrino emission from the 
"hot spot", which is marked by the point Pi 
(OPi = R) in Figure 1. At the point P (OP = 
R) nearby Pi, the ejection of matter is due to the 
local (lower) isotropic neutrino emission around 
Po. The matter suddenly sees a substantially 
larger number of neutrinos when passing through 
the point P2. Note that neutrino emission at an 
arbitrary point on the neutrino sphere (e.g., Po or 
Pi) is assumed to be isotropic in all directions (i.e., 
the local neutrino flux is not radial) as in Otsuki 
et al. (2000). This sudden increase of the neutrino 
number density at P2 is approximated by a jump 
of the neutrino luminosity from the original value 
L v = 1 x 10 51 ergss -1 for R < r < R 2 to the ef- 
fective luminosity L V 2 for r > R2, where r is the 
distance from the center O and R2 = OP2. 

The wind models considered in this study are 
listed in the first column of Table 1 , where R2 (sec- 
ond column) and L u2 (third column) are taken to 
be 10, 12, 15, and 20 km, and 1.0, 2.0, 3.0, 4.0, 
and 5.0 in units of 10 51 ergss -1 . The resulting net 
heating rate (q; top panel) for models Al, A5, B5, 
C5, and D5 and each heating/cooling rate (bot- 
tom panel) for models A5 and B5 are shown in 
Figure 2, as functions of r. Note that A1-A5 are 
isotropic wind models (i.e., R2 = R). 

For isotropic winds (A1-A5 in Table 1), the 
higher L V 2 (= L„ in these cases) results in shorter 
dynamic timescale T dyn (= \p/(dp/dt)\ T =o.5Mcv) 
and lower asymptotic (i.e., maximum) entropy s. 
This shows that the increased q (see Al and A5 
in Fig. 2, top panel) is consumed to drive more 
matter (i.e., higher M as can be seen in Table 1) 
from the neutron star surface with faster velocity, 
rather than to increase entropy. In contrast, for 
anisotropic models, an increase of L„2 (for r > R2) 
is quite efficient both to increase entropy and to 
reduce dynamic timescale (Table 1). The reason 
is that the matter has been already lifted with 
low L v (= 1 x 10 51 ergss _1 < L^) and thus with 
small M. Therefore, the density (and temper- 
ature) at arbitrary r is significantly small com- 
pared to the corresponding isotropic wind. This 
can be seen in the 5th (and 6th) column in Ta- 
ble 1 , which lists the density P13 (and temperature 
T13) at r = 13 km (see about one order difference 
in P13 for A5 and B5). 

For isotropic wind models, the five times greater 
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neutrino luminosity simply results in the increase 
of q with the same factor (Al and A5 in Fig. 2, 
top panel). This does not hold, however, for 
anisotropic wind models. For model B5, the maxi- 
mum q is as twice large as that for model A5 (with 
the same L v2 ), and more than 10 times larger than 
that for model Al (with the same L v ). This can 
be explained as follows. As shown in Figure 2 
{bottom panel), for isotropic winds {dashed lines; 
A5), the heating is mainly due to q u n and q ve , 
while q vv plays only a minor role. In contrast, for 
anisotropic winds {solid lines; B5), the neutrino 
pair annihilation q vv plays a crucial role, whose 
peak (at r « 13 km) is a factor of seven higher 
than that in A5. This effect can be clearly seen in 
Figure 2 {top panel), in which the case without an 
increase of q vv (model B5a) and with an increase 
of q vv only (model B5b) are compared (see also 
Table 1). 

This is due to the difference of p and T de- 
pendences in these heating terms. For a fixed set 
of r, Y e , L v and neutrino mean energies, these 
heating rates are related to p and T such as 
QvN = constant, q ve oc T 4 p _1 , and q vlJ oc p^ 1 
(Qian & Woosley 1996; Otsuki et al. 2000). As 
a result, q v ^ in B5 (Fig. 2, bottom panel) closely 
follows that in A5 for r > R 2 , which is indepen- 
dent of p and T. As can be seen in Table 1, a 
reduction in p owing to low L v in B5 (compared 
to that in A5) is accompanied with a reduction in 
T. As a consequence, q ve in B5 is lower than that 
in A5 even for r > R 2 . However, q vv is not depen- 
dent on T but is inversely proportional to p (i.e., 
proportional to the number of neutrinos per vol- 
ume), which becomes significantly high owing to 
the decreasing p for r > R 2 (Table 1). Note that 
the cooling terms {q e N oc T 6 and q ee oc T 9 /? -1 ) 
quickly decay with increasing r and have negligi- 
ble effects by the anisotropy (Fig. 2). 

As can be seen in the above numerical exper- 
iments, the strong anisotropy in neutrino emis- 
sion can be an additional energy source pointed 
out by Qian & Woosley (1996). However, this 
mechanism may work only for r < 1.5 R, which is 
rather closer to the neutrino sphere than the sug- 
gested range (1.5 < r/R < 3) by Qian & Woosley 
(1996). In the current study, the effects of increas- 
ing entropy and accelerating wind are prominent 
for the wind closer to the hot spot, in particu- 
lar for R w 12 km (Table 1), at which q maxi- 



mizes (Fig. 2). The effect of anisotropic neutrino 
emission becomes less important for a more dis- 
tant wind (e.g., R 2 — 20km), where the neutrino 
heating has mostly ceased (Fig. 2). For a fixed 
R 2 , the effect is more significant for higher L v2 as 
can be seen in Table 1. For model B5, the entropy 
is about 50% higher (180 A^ fc) and the dynamic 
timescale is about a factor of ten shorter (1.65 ms) 
than those in the isotropic model Al with the same 
L v {s = 117 NaIz and Td yn = 14.1ms). 

Note that the purely parametric examinations 
explored in this section should be regarded as only 
qualitative ones. For instance, the "hot spot" is 
not necessary a point as illustrated in Figure 1. 
It is conceivable that the area with strong neu- 
trino emission has some distribution on the neu- 
trino sphere. Moreover, the configuration should 
become multi-dimensional soon after the wind ma- 
terial passes the point P2 in Figure 1, which is 
treated within the framework of a spherical wind 
model in the current study. More realistically, the 
wind matter deviates from the radial to the di- 
rection of OP2 in Figure 1. This may moderate 
the acceleration of wind and increase the heating 
duration. Hence, the current results may over- 
estimate the reduction of dynamic timescale and 
underestimate the increase of entropy. It is diffi- 
cult to estimate the net effect to the nucleosynthe- 
sis from these modifications. Obviously, a multi- 
dimensional approach will be needed to quantita- 
tively estimate the effects of the anisotropy. 

3. Nucleosynthesis in Winds 

Adopting the wind trajectories discussed in § 2 
for the physical conditions, the nucleosynthetic 
yields are obtained by solving an extensive nuclear 
reaction network. The network consists of 6300 
species between the proton and neutron drip lines 
(for more detail, see Wanajo 2006). Neutrino- 
induced reactions and nuclear fission are not con- 
sidered in the current study. Each calculation 
is initiated when the temperature decreases to 
T 9 = 9 (where T 9 = T/10 9 K). The initial com- 
positions are given by X n = 1 — Y e and X p = Y e , 
respectively, where X n and X p are the mass frac- 
tions of neutrons and protons, and Y e is the initial 
electron fraction (number of proton per nucleon) 
at Tg = 9. In this study, Y e is taken to be 0.4, ac- 
cording to the core-collapse simulation in Woosley 
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et al. (1994, at L v « 1 x 10 51 crgss" 1 ). As in 
Wanajo et al. (2002), the temperature and density 
are set to be constant when Tg decreases to 1.0, 
in order to mimic the effect of the slower outgoing 
ejecta behind the shock. 

The nucleosynthesis results for models B2-B5, 
C2-C5, and D2-D5 (Table 1) are shown in Fig- 
ure 3, as a function of atomic mass number. For 
anisotropic wind models with R 2 = 12 km (B2- 
B5), the effect of anisotropic neutrino emission is 
evident. A factor of three or four increase in L u2 
(B3 and B4 in Table 1) leads to s w 150-160 N A k 
and Tdyn ~ 3 — 4 ms, resulting in the r-process nu- 
cleosynthesis (Fig. 3). For model B5, the high en- 
tropy (= 180 N a k) and short dynamic timescale 
(= 1.65 ms) of the wind drive the nuclear matter 
to the actinide region. The neutron-to-seed abun- 
dance ratio at the beginning of the r-process, de- 
fined as Tg = 2.5, is Y n /Yh = 176 and the fi- 
nal averaged mass number of heavy nuclei with 
Z > 2 is (A h ) = 230 (Table 1). For the models 
with R 2 = 15 km (C2-C5), the r-process still takes 
place when L v2 is four or five times higher than 
L v (models C4 and C5). For R 2 = 20 km (D2- 
D5), the effect of anisotropic neutrino emission is 
not important and the nucleosynthesis results are 
not significantly different from the isotropic cases 
(A1-A5). 

4. Implications 

In this Letter, the effects of anisotropy in neu- 
trino emission for the r-process nucleosynthesis 
in proto-neutron-star winds were examined, using 
the spherically symmetric, steady outflow model of 
neutrino-driven winds. It was shown that strong 
anisotropy, if it exists, can be an additional en- 
ergy source (Qian & Woosley 1996) to heat the 
wind material. A factor of four or five enhance- 
ment in effective neutrino luminosity results in 
the significant increase of entropy and shorten- 
ing of dynamic timescale of outgoing neutrino- 
heated ejecta. This is mainly due to the boosted 
neutrino heating from annihilation of ncutrino- 
antineutrino pairs into electron-positron pairs as 
a result of anisotropic neutrino emission. This 
provides the physical condition suitable for the 
robust r-process, producing the third abundance 
peak (A = 195) and beyond. 

It is conceivable that asymmetric neutrino 



emission can be associated with the anisotropic 
matter distribution near the neutrino sphere. As 
an example, Kotake et al. (2003) suggested that 
the non-spherical neutrino sphere owing to rapid 
rotation leads to anisotropic neutrino heating with 
the pole-to-equator ratio of a few to more than 
10. This may result in strong contrast in neutrino 
emission on the neutrino sphere, which forms an 
effective "hot spot" around the rotational axis. 
A recent work with more sophisticated neutrino- 
transport scheme by Walder et al. (2005) showed, 
however, that the pole-to-equator flux ratio is at 
most a factor of two, even for a rather rapidly 
rotating core. This is a consequence that the ra- 
diation field is smoothened by the many neutrino 
sources above the neutrino sphere (e.g., convec- 
tive bubbles) at the early phase (< 1 s after core 
bounce). Nevertheless, all the convective bubbles 
are evacuated during the late wind phase (~ 10 s) 
and a strong contrast of neutrino flux might form 
on the neutrino sphere for a rapidly rotating core. 

Another possibility of anisotropic neutrino 
emission might be due a global fluid instabilities 
of neutrino-heated matter as observed in multi- 
dimensional hydrodynamic simulations. Recent 
works have shown that hydrodynamic instabilities 
can lead to low-mode (/ = 1 in terms of an expan- 
sion in spherical harmonics of order I) oscillation 
of the convective fluid flow in the neutrino-heated 
layer behind the shock (e.g., Scheck et al. 2006; 
Buras et al. 2006b; Burrows et al. 2006). The 
presence of such a low convective mode results in 
the pair of a single outflow and a narrow accretion 
flow that creates the "hot spot" on the neutron 
star surface. It should be noted, however, that 
the two-dimensional simulations by Scheck et al. 
(2006) showed that the anisotropy of the accretion 
luminosity owing to this flow appears to be only 
a few percent (at least during the early phase up 
to ~ 1 s after core bounce) . A future investiga- 
tion relevant for the wind phase (~ 1 — 10 s) will 
be needed to examine the degree of anisotropic 
neutrino emission from such an accretion flow. 

Given one of the above (or another unknown) 
mechanism works, a constraint for the r-process 
may be obtained from the condition that creates 
the "hot spot" owing to, e.g., rapid rotation or 
long lasting accretion flow. It is conceivable that 
only a limited fraction of supernovae create the 
"host spot" relevant for the current scenario (e.g., 
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rapid rotators or less-energetic supernovae that 
form the long lasting accretion flow) . This can be 
a reasonable explanation for that the spectroscopic 
analysis of extremely metal-poor stars and Galac- 
tic chemical evolution study imply only a limited 
fraction of core-collapse supernovae undergo the re- 
process nucleosynthesis (Ishimaru & Wanajo 1999; 
Ishimaru et al. 2004). 

The implications in this study must be tested 
by future multi-dimensional simulations of core- 
collapse supernovae for long duration (~ 10 s) with 
accurate neutrino transport. Systematic calcula- 
tions of nucleosynthesis with such hydrodynamic 
trajectories will be also needed to investigate the 
contribution to the Galactic chemical evolution of 
r-process nuclei. 
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Table 1 

Model Parameters and Results 





i?2 


Lu2 (1C 


I 51 
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P13 (10 7 


Tia (10 10 


s 


Tdyn 






Model 


(km) 


ergs s~ 




MqS- 1 ) 


g cm -3 ) 
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(N A k) 


(ms) 


Y n /Y h 


(A h ) 


Al 


10 


1.0 




3.86 


3.49 


2.06 


117 


14.1 


6.01 


107 


A2 


10 


2.0 




13.3 


5.98 


2.35 


103 


7.20 


8.16 


109 


A3 


10 


3.0 




27.4 


8.10 


2.54 


95.5 


4.99 


6.02 


107 


A4 


10 


4.0 




46.0 


9.75 


2.66 


90.7 


3.90 


6.02 


107 


A5 


10 


5.0 




68.7 


11.9 


2.80 


87.0 


3.26 


10.3 


111 


B2 


12 


2.0 




4.88 


2.68 


1.88 


131 


7.45 


18.0 


118 


B3 


12 


3.0 




5.84 


2.08 


1.73 


145 


4.44 


38.9 


135 


B4 


12 


4.0 




6.84 


1.65 


1.61 


161 


2.72 


80.7 


170 


B5 


12 


5.0 




7.93 


1.15 


1.46 


180 


1.65 
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B5a a 


12 


5.0 




5.84 


2.07 


1.74 


147 


4.37 


39.8 
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B5b b 


12 
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2.61 
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174 
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15 


2.0 




4.15 
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127 


9.76 


12.1 
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C3 


15 


3.0 




4.41 


3.25 


1.94 


136 


7.08 


20.9 


120 


C4 


15 


4.0 




4.68 


3.25 


1.90 


147 


5.15 


34.7 


131 


C5 


15 


5.0 




4.97 


3.10 


1.84 


159 


3.69 


58.3 


153 


D2 


20 


2.0 




3.93 


3.45 


2.04 


122 


11.9 


8.43 


109 


D3 


20 


3.0 




4.00 


3.53 


2.04 


127 


10.2 


11.3 


112 


D4 


20 


4.0 




4.08 


3.50 


2.02 


132 


8.69 


14.9 


115 


D5 


20 


5.0 




4.15 


3.48 


2.01 


137 


7.43 


19.3 


119 



a Same as B5, but without enhancement of q vv for r > i?2- 
b Same as B5, but with enhancement of only q vv for r > J?2- 
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Fig. 1 . — Illustration of asymmetric neutrino emis- 
sion. O is the center of the neutron star. Strong 
neutrino emission from the "hot spot" near the 
point Pi on the neutrino sphere is assumed, oth- 
erwise being isotropic. A wind blowing from a 
nearby point Po with the (weaker) isotropic neu- 
trino emission (L v ) suddenly see a larger number 
of neutrinos {L V 2) when passing P2. 




10 15 20 25 

radius (km) 

Fig. 2. — Top: Net neutrino heating rates for wind 
models Al, A5, B5, B5a, B5b, C5, and D5 listed 
in Table 1, as a function of r. Jump of the heat- 
ing rate at r = i? 2 is due to the sudden increase 
of effective neutrino luminosity from L v to L I/2 - 
Bottom: Heating (q V N, q_ve, and q vv ) and cooling 
(q e N and q ee ) rates as functions of r. Dashed and 
solid lines are for wind models A5 and B5 (listed 
in Table 1), respectively. 
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mass number mass number mass number mass number 

Fig. 3. — Final abundances obtained by the nu- 
cleosynthesis calculations for wind models listed in 
Table 1 (except for A1-A5) as a function of atomic 
mass number. 
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